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Abstract: A certain sector of the matrix model for M-theory on a plane wave background 
has recently been shown to produce the transverse five-brane. We consider this theory at 
finite temperature. We find that, at a critical temperature it has a Gross- Witten phase 
transition which corresponds to deconfinement of the matrix model gauge theory. We 
interpret the phase transition as the Hagedorn transition of M-theory and of type II string 
theory in the five-brane background. We also show that there is no Hagedorn behaviour 
in the transverse membrane background case. 
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1. Introduction and Summary 

The matrix theory conjecture states that the large N limit of a certain matrix quantum me- 
chanics is an exact description of light-cone quantized M-theory in flat eleven-dimensional 
space-time (lj (for a review see Q). Many of the objects of M-theory, such as the su- 
pergraviton, supermembrane and M five-brane emerge in a natural way from the matrix 
model and some of their low-energy interactions have been computed and matched with 
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supergravity. However, the study of this seemingly simple model has been plagued by its 
nonlinearity and the existence of flat directions in the potential leading to a continuous 
spectrum |J. 

Recently a matrix model of M-theory in an 11-dimensional plane- wave geometry has 
been studied P]-JlO[, It has the great advantage that one can formulate an accurate per- 
turbative expansion of the model. Properties of states and symmetries have been analyzed 
in this limit H-R- It has also been realized that some of the states are protected by 
supersymmetry. 

The plane-wave background has a single dimensional parameter fj, which appears in 
the metric and the constant four-form field strength, 

ds 2 = -2dx + dx~ " ((f) 2 {x a f + (f ) 2 (x a ') 2 ) (dx + ) 2 + {dx a f + (dx a ') 2 
F+123 = il (1-1) 

This is a 1-parameter family of backgrounds which goes to Minkowski space when fi — > 0. 
The matrix model for M-theory on this background has action [|J] 

S = RjdtTi^ (D t X'f + ^l>'D,i> + i,V [jP,</>] + i [X'.X'f 

~\ (if <*"> 2 " \ (if ~ 'i^' ~ ^-™} (") 

We have set the M-theory Planck length to one, £ p = 1, and the conventions for indices are 
... = 1, . . . , 9 and a, b, . . . = 1, 2, 3, a', b', . . . = 4, . . . , 9. All degrees of freedom are 
N x N Hermitian matrices. The covariant derivative is 

D t X = j t X + i[A,X] (1.3) 

This matrix model is conjectured to be a discrete light-cone quantization (DLCQ) of M- 
theory where R is the compactification radius of the null coordinate x~ and there are 
./V units of light-cone momentum, p + = N/R. It is a 1-parameter deformation of the 
BFSS matrix model Q which described M-theory on Minkowski space. Just as the BFSS 
model could also be regarded as the effective action describing the low energy dynamics 
of a collection of N DO-branes of type IIA string theory on Minkowski space, with the 



appropriate re-interpretation of the parameters, (1.2) could also describe DO-branes in a 
background of type IIA string theory. 

The matrix model can be analyzed in perturbation theory. As was noted in ref. Q, 
one interesting feature of the semiclassical expansion is the existence of a large number of 
different vacua. The potential contains a term of the form 

,JJLx a + if abc x b : 

Classical vacua are 



V ~ Tr ^X a + ie abc X b X c ) (1.4) 
V 3R 



x a (t) = (i.5) 
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where J a are matrices which satisfy the SU (2) Lie algebra, 



J a ,J b =ie abc J c (1.6) 

and could be either irreducible or reducible representations of 577(2). The representa- 
tions of SU(2) are interpreted as fuzzy spheres which are spherically symmetric states of 
membranes j|] (see also || ||). They approach a classical membrane when the spin of an 
irreducible representation is large, the state with a single membrane being described by a 
single TV-dimensional irreducible representation of SU(2) with spin j = ^p-. 

Recently, an alternative interpretation of these solutions as transverse five-branes has 
been given [1C]. In principle, any of the classical vacua corresponds to a five-brane state. 



A single classical five-brane, with its five spatial dimensions forming a five-sphere and its 
time direction lying along the light-cone, is believed to correspond to the vacuum where 
the representation of SU{2) is highly reducible, consisting of N repetitions of the singlet 
J 1 = representation. A detailed comparison of the fluctuation spectrum of the matrix 
model about these states with the fluctuations of the classical five-brane on a plane-wave 



background was given in ref. [ 10 1 and good agreement was found. A stack of k coincident 
five-branes is thought to correspond to the state with n = N/k repetitions of the k- 
dimensional representation of SU(2). If the representations have differing dimensionalities, 
with maximum dimension k, this still corresponds to k spherical five-brane states where 
the branes have different radii. 

The main difference between a five-brane state and a membrane state is that, for the 
five-brane, as N is taken to infinity, the number of repetitions of irreducible representations 
of SU(2) is taken to infinity with their dimensionalities held fixed, whereas for the mem- 
brane, the number or repetitions of an irreducible representation (number of membranes) 
is held fixed and the dimensionalities are taken to infinity. 

The distinction between these two cases that will be important to us in the following 
occurs in the way that the gauge symmetry of (|1.2| ) is realized. Vacua with non-trivial 
representations of SU (2) break the gauge symmetry, realizing it in a Higgs phase. There is 
a residual gauge symmetry which interchanges representations of the same dimensionality. 
If there are representations of dimension k, the gauge symmetry breaking pattern is 



n (It) <-) 



U(N) 

Here, because all degrees of freedom transform in the adjoint representation, the true gauge 
group is the group modulo its center. For a membrane state, as we take the limit N — ► oo, 
the rank of the residual gauge group remains finite, whereas for a five-brane, it becomes 
infinite. 

The coupling constant which governs the size of corrections to perturbation theory 
depends on the classical vacuum about which one is expanding. For membrane vacua, the 
effective coupling constant is 

9cs= {W (L8) 
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For the five-brane, where the multiplicity of representations is n, the effective coupling is 
the 't Hooft coupling 

g 2 eS n = A = n (1.9) 

The difference comes from the residual U(n) symmetry of the five-brane solution and the 
fact that index loops in perturbation theory contribute factors of n which can only be 
controlled in the 't Hooft limit, 

n —* oo with A fixed (1-10) 
The decoupling limit, which isolates the theory living on the five-brane, is the same as 



the 't Hooft limit in ( 1,10 ). Intuitively, this is the weakest coupling limit that can be taken 
that retains any of the interactions in the matrix model. If it were taken for a membrane 
state, the matrix model would be non-interacting. It is only the five-brane states which 
retain interactions. The 't Hooft coupling is the radius of the five-brane squared in units 
of the string scale 

(1.11) 

a 

In this work we will be concerned with perturbative computations in five-brane back- 
grounds in the 't Hooft limit. Our perturbative analysis is limited to small A, the situation 
where the five-brane is highly curved. In the end this will prevent us from making a direct 
comparison of our results with the properties of the usual NS five-brane of string theory 
which is large and flat and which we believe would be obtained by the infinite radius, and 
therefore infinite 't Hooft coupling limit, in these sectors of the matrix model. 

In this paper we shall analyze the thermodynamics of this matrix model. The main 
observation will be that the five-brane states of the matrix model have a first order phase 
transition at some critical temperature. This phase transition is seen in the 't Hooft 
limit. Our analysis is further limited to the regime where the 't Hooft coupling A is small. 
However, the existence of the phase transition appears robust and is expected to occur over 
a range of A. 

We will show this in two different ways. The first way counts the degeneracy of states 
with high energies in the infinite N limit and establishes that there is a Hagedorn density 
of states. From the exponential behavior we can find the Hagedorn temperature. 

Then, secondly, we use matrix model techniques to show that there is a Gross- Witten 
type of large N phase transition in the perturbative limit of the model. In the weak 
coupling limit, the transition temperature coincides with the Hagedorn temperature that 
we extracted from the asymptotic density of states. The advantage of the matrix model 
approach is that it is at least in principle possible to study the phase above the Hagedorn 
transition. The transition is similar to deconfinement transitions which are expected to 
occur in confining gauge theories at high temperature. 

Note Added 

While this research was conducted, an extensive study of the Hagedorn transition in the 
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case of weakly coupled SU(N) gauge theories on compact spatial manifolds was carried 
out in an interesting work |36|]. This work has a partial technical overlap with ours. 



2. Hagedorn Behaviour of Matrix Models 

In this section, we will compute the Hagedorn temperature of matrix models by counting 
their gauge invariant states, or, equivalently, the numbers of gauge invariant operators. We 
will show that, when the dimension N of the matrices is infinite, they have a Hagedorn 
density of states. Our analysis is very similar to that in ref. [11, [l^, 13]. 

2.1 A Toy Model 

We will begin with a toy model which captures the salient features of the physics involved. 
Consider a set of d bosonic N x N Hermitian matrices, Xj(t), where j = 1, . . . , d and t is 
the time. These are coupled to an N x N "gauge field" A(t). The action is 

S = j ^E^((^i+^»^]) 2 -^|) (2-1) 

If it were not for the gauge field appearing in the covariant derivative this would be a 
theory of d independent matrix oscillators. Introducing the gauge field makes the model 
invariant under the gauge transformation 

Xj(t) -» U(t)Xj(t)U\t) 
A(t) -» U{t)A{t)U ] {t) - iU(t)U\t) (2.2) 

The equation of motion which is obtained by varying A(t) has no time derivatives and is 
an equation of constraint 

d 

g = Y f [X j (t),IL j (t)]~0 (2.3) 

3=1 

where 

K 3 (t) = j t X j {t)+i[A{t),X J {t)] (2.4) 
is the canonical momentum, with Poisson bracket 

{xf(t),H% d (t)} = 5 jk S ad S bc (2.5) 



The left-hand-side of the constraint equation (^j) is the generator of an infinitesimal time- 
independent gauge transformation, 

SXj = {Tr (AG) ,Xj} = - [A, X,] (2.6) 

When we quantize, this toy model is simply described by d decoupled matrix oscillators 
with the constraint on the physical states. 1 



1 Obviously, in a deconfined phase this no longer holds for the vacuum as the center symmetry is broken. 
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Claim. In the limit N —* oo this toy model has a Hagedorn density of states and Hagedorn 
temperature 

T H =u/\nd (2.7) 

The case d = 1 should be interpreted as Th = oo. We will also use the parameter (3 = 1/T 
for the inverse temperature and the notation f3jj = 1/Tjj = Ind/u. 

To substantiate this claim, we will count the number of physical states which have a given 
asymptotically large energy and extract f3jj from the exponential growth (if there is any). 
The Hamiltonian which follows from canonical quantization of the action (|2.1] ) is 



H = ^T±Tr{n* + uj 2 X]) (2.8) 

where the momentum obeys the canonical commutation relation 

Xf ,Wjf] =iS jk 5 ad 5 bc (2.9) 



It is easy to see that the constraint (2.3) commutes with the Hamiltonian. 

We shall take the strategy of quantizing the model of free oscillators with Hamiltonian 
Q2.8D and commutators (|2.9| ) . From the resulting quantum states we then isolate a subspace 
of physical states which are annihilated by the constraint operator in equation ( |2.3[ ). 

We define creation and annihilation operators 

' (IL - iujXj) , at = — == (Uj + iujXj) (2.10) 



so that 

ab J cd ] _ x . x ad x bc 



a j > a k 



5 jk 6 aa 5 0c (2.11) 
and, after dropping the ground state energy, the Hamiltonian is 

d 

H = ^ w Tr (a)a 3 ) (2.12) 
i=i 

The vacuum state is annihilated by all of the annihilation operators 

af\0)=0 Va,6 = 1, ...,N;j = l,...,d (2.13) 
A basis of quantum states is found by operating creation operators on the vacuum, 

|0) , af |0) , afal cd \0)... (2.14) 



These basis vectors are eigenstates of the Hamiltonian. 

From these states, we must choose a subspace which is annihilated by the constraint 
operator in (|2.3j). That operator generates infinitesimal time independent gauge transfor- 
mations. Quantum states will be invariant if they are invariant under the finite transform 

aj -► UayU* , a] -► Ua)u ] , with UU ] = 1 (2.15) 
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Such invariant states are generally obtained from the general basis states by taking traces 
over the matrix indices, 2 

"h in\ tv /^"LtA m\ Tv ( r,^\ tv ( /it 



|0) , Tr^ T j|0) , Tr[a)a\)\U) , Tr (a) J Tr (a\) |0) . . . (2.16) 
The operators which create gauge invariant states are thus of the form 



Tr ( at . . . at 



11 Jn t 

and can be thought of as "words" which are made from the "letters" a\. The length of a 
word is the number of letters which it contains. There are d different letters at our disposal. 
Inside a word, the ordering of the letters is important, since they are matrices and do not 
commute with each other. However, because of the cyclicity of the trace, words which are 
related by cyclic permutations of the letters are equivalent. 

Note that for finite iV the words are not independent. Generally, a very long word 
can be expressed in terms of linear combinations of sentences made from shorter words: in 
the U(N) theory, the operators Tr(a^) fc for k > N can be expressed in terms of those with 
< k < N. However, the long words do become independent when the limit N — ► oo is 
taken. Here, we shall always assume that the large N limit is to be taken. For this reason, 
we shall treat all distinct words as independent operators. 

Finally, we can create a state by applying several "words" to the vacuum. Such multi- 
trace operators, 

Tr ( a h • • • a l) Tr (°L • • • °L) • • • Tr ( a l ■ ■ ■ a l p 
will be called "sentences" . The length of a sentence is the sum of the lengths of its con- 
stituent words, or the total number of letters in the sentence. As each creation operator, 
or letter, creates an excitation of energy uj, the energy is the length of the sentence in units 
of to. 

Obviously, words or sentences should be considered equal if they produce the same 
quantum state. As operators, words commute with each other. Therefore, the order of the 
words in a sentence does not make any difference. 

In order to find the Hagedorn temperature we should find the density of states at 
asymptotically large energy. If the theory has a Hagedorn spectrum, the increase in the 
density of states will be exponential, 

p(E) ~ e fe£ (2.17) 
This formula can be corrected by prefactors; we shall discuss them later, at subsection 



2.2 A. We will concentrate for the time being on extracting the coefficient in the exponent 
which we shall identify with the inverse of the Hagedorn temperature. 

To find the density of states, we are essentially looking for the number s n of sentences 
of length n. In order to compute this number, we will need to consider also the number 
Wk of words of length k. It is convenient to define 

so = w = 1 (2.18) 



2 If N were finite there would also be the possibility of other invariants such as the determinant of aj. 
The energy of the state created by this operator is Nu. 
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As is common in combinatorial problems, it is helpful to define the generating functions of 
those numbers, 

oo 

W(x) = ^2w k x k (2.19) 

fc=0 

oo 

S(x) = J2 s nX n (2.20) 

n=0 

2.1.1 The d = 1 Case 

For the case d = 1 of an alphabet with a single letter a) = a{, there is obviously only one 
word Tr(a^) fc of each length k, or w k = 1. A sentence of length n is nothing but a collection 
of such words, where their order is immaterial. Hence, the number of sentences of length 
n is just the number of partitions of n into natural numbers, s n = p(n). It is well known 
from the days of Euler that 



oo 



P(x) = ^2p(n)x n = (1 + x + x 2 + . . .)(1 + x 2 + x 4 + . . .)(1 + x 3 + x 6 + ...)•• • 

n=0 

oo 

= H(l-x i r 1 (2-21) 

and from the celebrated work of Hardy and Ramanujan that 

p(n) ~ e c ^ (2.22) 

with c = ~ 2.5651. 3 

The number of sentences of length n is therefore sub-exponential, and the partition 
function converges for all finite temperatures; consequently we conclude that (3h = and 
in this case Tjj = oo as claimed. 

2.1.2 The d > 1 Case 

An upper bound on the number of words of length k is 

Wk < d k 

as there are d independent choices for each letter. The cyclic property of the trace makes 
some of those words equivalent. Obviously, each word can have at most k equivalent 
representations; therefore, a lower bound on the number of inequivalent words is 

w k >\d k 
k 



3 



In fact, the result is much stronger. In particular 



p(„) = * ( 1 - J-) e^ + o( -e^ 2 



where 

1 

n = n 

24 
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We claim that the latter bound is in fact a good approximation, 

w k « -d K 

To substantiate this claim, we will utilize the famous Polya theorem. 

Theorem (Polya). Let G be a group acting on a set X . Then the number of colourings 
of the elements of X with d colours, distinct under the action of G, is 

where X(g,X) is the number of cycles in the action of g on X. 

In our case, X = Z/% is the string of letters in a word, and the group of rotations is also 
G = Zfc, acting on X by addition. The colours are the different kinds of letters, so their 
number is indeed d. The number of colourings is just the number of inequivalent words. 

Let us begin with the case when k = p is a prime number. Then, the identity has 
\X\ = p orbits, and each other rotation has only one; therefore, 4 

w p = X - (d p + {p- l)d) (2.24) 

This is indeed very close, for large p, to |<i p . For the case where k is not a prime, a similar 
formula can be written in terms of Euler's totient function, but it suffices to note that apart 
from the identity, again having \X\ = k orbits, all other rotations can have at most k/2 
orbits (in fact, there can be at most one such orbit, for g = k/2 when k is even). Therefore 

w k < i (d k + (k - l 

which is again close to \d k . 

Accordingly, we shall henceforth use the approximation 

w k w -d k 

for k > 1. This results in 

W(x) = w k x k w 1 + -^d kxk = 1 - lo g(! - dx ) ( 2 - 25 ) 
k=0 k=l 

A sentence can be composed of only one word; The contribution to s n is then w n , 
and to S(x) is W(x). It can be composed of two words; the contribution to s n is then 
naively the convolution (w * w) n = J2m=o w mW n - m , and to S(x) is W(x) 2 , but we should 
remember that the order of the words is immaterial. As for large n the sentences composed 
of two identical words are negligible (the generating function of those is W(x 2 )), we should 



incidentally, the fact that w p is integer provides a simple proof of Fermat's little theorem, <P 
(mod p). 
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approximately divide by two. In general we approximate that the number of words is 
so large that generically no two will be identical in a sentence; if there are a words in a 
sentence, we should divide the naive approximation by a\. Consequently, 

oo 1 ^ oo 

S{x) « -j W ( x ) a = eW(x) ~ e • Y^dx" = e '^ ^ (2 ' 26) 

a=0 ' n=0 

SO 

s n w e • d n ~ e ( logd ) n (2.27) 

and the Hagedorn temperature, in units of u;, is indeed as claimed. 
Note that if 

x = e-P" (2.28) 

then 

oo 

Z(/3) = S(x) = Y,s n e- pm} (2.29) 

ra=0 



is exactly the partition function, so in light of ( |2.26| ) we get 

z( «~ e 'T3^ (2 - 30) 

The partition function clearly diverges at the inverse Hagedorn temperature (5h such that 

d . e -P H ^ = i (2.3i) 

which gives indeed (3h = 

In fact, a more careful use of Polya's theorem can give an exact expression [11] for 
Z(P) or S(x), but this is irrelevant for our purposes. 

2.2 Extensions 

2.2.1 The Traceless Case 

The matrices we have encountered thus far were assumed to be Hermitian. They can 
equivalently be characterized as matrices in the adjoint representation of the U(N) algebra. 
It is natural to explore the variant where the matrices are in the adjoint of SU(N) (we will 
denote the various magnitudes for the SU(N) version with a tilde). We will now show that 
our results are essentially unchanged if we have the gauge group SU(N) instead of U(N). 

A matrix in the adjoint representation of SU(N) is Hermitian, but with the extra 
constraint of being traceless. In particular Tra^ = and so there are no single-letter 
words: 

wt = (2.32) 

In the d = 1 case, therefore, the partitions of n should not contain ones. As any 
partition without ones of any m < n can be uniquely completed, with the addition of 
n — m ones, to a partition of n, and vice versa, obviously 

n 



~ s m (2-33) 



m=0 
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and so 

d 1 . 

s n = s n - s n - 1 pa — s n ~ —j=s n (2.34) 
on Jn 



We see that the exponential behaviour, and hence the Hagedorn temperature, remains as 
in the U(N) case, and the change is only in the irrelevant polynomial prefactor. All this 
can also be seen from the fact that the relevant partition function is 

oo 

P{x) = Y[(l - xT 1 = (1 - x)P{x) (2.35) 
i=2 

In the d > 1 case, ( p. 32 ) implies 

W(x) = W(x) - dx (2.36) 

and 

S{x) pa e^(*) = e ~ dx e w ^ pa e~ dx S{x) 

pa e (1 - dx + -d 2 x 2 - -d 3 x 3 ± . . .) (1 + dx + d 2 x 2 + d 3 x 3 + . . .) (2.37) 

so the coefficient of x n is now 

Sft „ ed »(l-l + i-I + . .. + (-!)»!) 

2 3! n! 

but the last factor is approximately e , so 

s n pa d n (2.38) 

We get the same Hagedorn temperature as in the U (N) case, the only change being in an 
irrelevant constant factor. 

2.2.2 Fermions 

Fermionic matrix degrees of freedom behave quite similarly to bosonic ones of the same 
mass. We can form the fermionic creation operators b\. Each entry of such a matrix is 
a Grassmannian number which squares to zero; however, there are so many of those, (i.e. 
iV 2 ) that when forming words of finite length, none vanishes, and more generally, the words 
remain independent. However, each word squares to zero, so we can not use the same word 
more than once in a sentence. 

For the d = 1 case, this means that the partitions of n should be to distinct numbers. 
The number of such partitions, q(n), behaves also as 

q(n) ~ e c '^ (2.39) 

when now d = \J\k. The corresponding generating function for the sentences is easily 
seen to be 5 



° Incidentally, the number q(n) of partitions of n into distinct parts is equal to the the number of partitions 
into odd parts; this is an amusing exercise in generating functions, as formally 

oc oo 

fl(l + x i ) = n(l-^" 1 )" 1 - (2-40) 

i=l 3=1 
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S(x) = Q(x) = Q(n)x n = H(l + s*) (2.41) 

n=0 8=1 

For the d > 1 case, we should again take into account the restriction that no two equal 
words appear inside a sentence. However, we have already argued that such cases are rare 
and used the approximation that they are negligible, when we argued that the "a words in 
a sentence" term W(x) a should be divided by a! (so that S(x) is the exponent of W(x)). 
We conclude that the Hagedorn temperature of the d > 1 fermionic matrix model is equal 
to its bosonic counterpart. 

2.2.3 Different Masses 

When there are several matrices with different masses, there is obviously still a Hagedorn 
density of states, with the Hagedorn temperature between two bounds: the Hagedorn 
temperature of the same number of matrices but all with the same mass, which is either 
the minimal or maximal one. However, it is quite easy to generalize (2.31) to get an exact 
result. 

Claim. For a model with d > 1 matrices, bosonic or fermionic, of masses LOi > 0, i = 
1 . . . d, the inverse Hagedorn temperature Ph is given by the solution of 



= 1 (2.42) 



i=l 

To see that, we first pick a mass scale uj and write 

u>i = ujvi (2.43) 
The idea of the proof in the equal masses case carries over when we define 



F{x)=^jT l x Vi (2.44) 



i=l 

as a generalization of dx. We still have 



CO . 

W(x) ss 1 + j;F(x) k = 1 - log(l - F(x)) (2.45) 

and 



k 

k=l 



S(x) * e w ^ „ e . _1_ (2.46) 

Where now s n should be interpreted as the number of states with energy noj and not as 
the number of sentences of length n; indeed, n doesn't even have to be integer. 
Clearly, S(x) diverges at 

x H = e~ tow (2.47) 

where 

F(x H ) = 1 (2.48) 
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Still using (12.28|) , we get 



Z{(5) 



(2.49) 



1 - Eti e-** 

which diverges at the claimed inverse Hagedorn temperature. As F is a monotone decreas- 
ing function of f3, where F = d > 1 for (3 = and i 7 = < 1 for (5 = oo, there is exactly 
one root of equation ( |2.42| ). 

If uJi = for some i, then there is an infinite number of zero modes, the Hagedorn 
temperature is formally zero, and we see that the relation (2.42) continues to hold formally. 

It is also easy to see that this reasoning, properly interpreted, does not rely upon the 
number of matrices being finite. 

2.2.4 Absence of Prefactor 

We would also like to refine our analysis and find any energy power prefactor multiplying 
the leading exponential in the density of states. Stated otherwise, assuming that the density 
of states behaves as 

p(E) ~ C 'E a e^ E (2.50) 

we wish to find a. 

Claim. There is no power-like prefactor in matrix models, i.e. a = 0. 

In order to see this, assume ( |2.50| ), and compute the partition function, 

Z{fi) = [ p(E) e-? E dE = T C E a e-^~M E dE = 

J Jo (P-M a - 

or equivalents by ( ^28| , ^29| , ^4^ ) , 



+i 



(2.51) 



S(x) 



CT(a + l)uj a+1 



It is easy to see that 



(logixn/x))^ 1 
S(x)S"(x) a + 2-log(x H /x) 



(2.52) 



(2.53) 



S'(x) 2 a + 1 

so that q can be extracted from the behaviour of the latter expression at the divergence, 



S(x)S"(x) 



S'(xY 



x=xh 



a + 2 
a + 1 



On the other hand, it is immediate from ( p.46| ) that 

S(x) S"(x) 2F'{x) 2 - (F{x) - l)F"{x) 



S'{x) 



F'{xf 



and therefore, from (2.48), 



S{x) S"{x) 



(2.54) 

(2.55) 
(2.56) 



X=X H 



S'(x) 2 

Equating ( ^.54 ) and ( 2.56 ) we get indeed that a = as claimed. 

Essentially, this behaviour is dictated by the fact that F(x), a monotonous function 
has a non zero derivative at x = xh- 
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2.3 Hagedorn Behaviour in Free Second Quantized Theories 

Let us study a free quantum theory of harmonic oscillators of frequencies u>i > 0. Their 
number may be finite or infinite, and they might be bosonic or fermionic. The correspond- 
ing Hamiltonian, after discarding the zero-energy, is H = ^^Wjajaj. The Fock space 
states of the second quantized theory (the "sentences" ) are generated directly by applying 
the creation operators (the "letters") on the vacuum: as the creation operators are not 
matrices, and no gauge invariance condition is imposed, there is no need for "words" in 
this context. 

Claim. The inverse Hagedorn temperature (3h in this model is given by the solution of 

^V &£ ^ = oo (2.57) 

i 

More accurately, it is the maximal value of (5 where the left hand side diverges. 

Note that this implies that the Hagedorn temperature for a system with a finite number 
of oscillators is always infinite, so there is no Hagedorn transition (for that, it is essential 
that there are no zero modes oJi = 0) . Even for systems with an infinite number of oscillators 
the Hagedorn temperature might be infinite. It is indeed infinite in the sensible cases where 
the model comes from a field theory in a finite volume. 

To understand the claim, we will begin with the example of d bosonic oscillators of 



equal frequency ui. As in the the corresponding matrix case of subsection [2.1. 4 s n counts 
the number of sentences of length n made out of d letters. The order of the letters does not 
matter, as the creation operators commute, and therefore s n = _7 1 ) (think of arranging 
n circles and d — 1 dividers in a line). As the density of states s n ~ const • n d ~ l behaves 
polynomially and not exponentially, there is no Hagedorn behaviour. 
This can also be easily seen from the generating function: 

W = t( n+ /_~ 1 i y= ah¥ = A (1 + x + x2 + x " + ■ ■ } (2 - 58) 

n=0 ^ ' ^ ' i=l 

The term x 3 in the i-th multiplicand represents j excitations of the i-th. oscillator. 

The generalization to different frequencies is immediate: using ( p. 28 , 2.43), a bosonic 



oscillator a) contributes a factor of ^ „ , while a fermionic one can be excited at most 
once, as (b^) 2 = 0, and so contributes a factor of 1 + x v . Accordingly, 

s( X )= n y^i- n (1+**) ^ 

Bosons Fermions 

Now, take the logarithm: 

\ogS(x)= -log(l-a;^)+ ^ log(l + x^) (2.60) 

Bosons Fermions 

If there is a subsequence of the going to zero, then log S(x) diverges for all x > 0, so 
by ( 2.28| ), [3h = 00; this is clearly true also from ( |2,57| ). Assume, on the other hand, that 
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u* = inf Vi > 0. Then for 1 > x > we have x Vi < x u * < 1, so there is some r\ < 1 such 
that (1 - r))x n < - log (1 - x n ) , log (1 + x Vi ) < (1 + rj)x n . Consequently, ( p.60| ) and the 
rig ht hand of ( ^57|) converge or diverge together. This concludes the proof of the claim. 

As an example, let us take a classical (non-relativistic) string tied in its two ends and 
free to oscillate in one transverse dimension. The string's frequencies are harmonics of 
some natural frequency, u>i = i ■ co for i = 1, 2, 3, . . .. The number of states of energy n ■ u> 
is again p(n), the number of partitions of n, which is sub exponential ( |2.22j ), so there is no 
Hagedorn behaviour. S(x) = P(x) is given by (2.21) that has a radius of convergence one 
around the origin. In this case, The right hand side of ( |2.57| ), ]T\ e _/3a;i = ^2, i x % = 
has indeed the same radius of convergence. If the string lives in D spacetime dimensions, 
there are D — 2 transverse directions, so S(x) = P(x) D ~ 2 , which still displays no Hagedorn 
behaviour. 

The lack of Hagedorn behaviour is generic in field theories. In the free case in a finite 
volume, the number of oscillators grows like the phase space, or as a power of the energy. 
The result ( |2.57| ) implies that we should look at expressions of the form i d x % ~ 

; as was shown, in a different context, in ( |2.58| ). The divergence is still at xh = 1, 
corresponding to (3h = or Tjj = oo. 



3. The BMN Matrix Model 



As mentioned in the introduction, this model arises from the study of a certain M-theory 
plane wave which is the Penrose limit of the AdS4 x S 7 or AdSy x 5* 4 solutions. An 
interacting matrix model arises from the DLCQ description of this background when there 
are units of momentum on the light-like circle [Q] . 

There are 11 — 2 = 9 coordinates apart from the light cone ones, which break into two 
groups: three coordinates which are equivalent under an SO(3) rotation, and the other 
six, equivalent under SO (6). Those coordinates give rise to nine bosonic matrices: three 
of mass ^oj for the appropriate oj, and six of mass gw. The model is supersymmetric and 
correspondingly has eight fermionic matrices, all of mass -jOJ. 

The various classical solutions of the matrix model are given by vanishing SO(6) (and 
fermionic) matrices, and by the SO(3) matrices being a dimension TV representation of 
SU(2). Taking the small 't Hooft limit makes the theory free around the classical solution. 



3.1 A Single Five-Brane 

A single spherical five-brane solution is represented by the classical solution with the SO(3) 
matrices vanishing (or being in the trivial representation) |uj. This vacuum manifestly 
respects the U(N) gauge symmetry of the model. The oscillators are then just given by 
the bosonic and fermionic U(N) matrices described above. 



According to our previous result (2.42), the inverse Hagedorn temperature is given by 
the solution of 

3e -/V3 + 8e~^/ 4 + 6e~^/ 6 = 1 (3.1) 
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or, if x = e ^/ 12 ; by the solution < x < 1 of 

(x + lf{x --) = 3f 4 + 8x 3 + 6x 2 - 1 = (3.2) 
3 

that is, by x = §, so /? H = - logx • ^ = « or T# « 0.0758533a;. 

3.2 Multiple Five-Branes 

A solution of k coincident spherical five-branes corresponds to the matrix model classical 
vacuum with the SO (3) matrices being n copies of the k dimensional representation, where 
n = N/k is taken to infinity, with k fixed. The SO(3) matrices having a vacuum expectation 
value results in the Higgsing of the gauge group to U(n). Each U(N) matrix is decomposed 
into k 2 different U(n) matrix oscillators, which arrange themselves into multiplets of SO(3) 
spins j, integer for the bosonic oscillators and half integer for the fermionic ones. As k 
grows, more and more oscillators are added, as the maximal j grows by one for each 
oscillator type. Essentially, each oscillator becomes a tower of oscillators, with the original 
k = 1 oscillator at the bottom, and two new towers begin at k = 2. As a consequence, 
the Hagedorn temperature is lowered. We will see, however, that the change is very small 
even in the large k limit, because the masses of the oscillators grow linearly in j (their 
degeneracy does also, of course, but this is much less important). 

The spectrum of oscillators in this case was computed in [||, g], and is given (in units 
of u) in table |], where we also give the labels given there to the different oscillator towers. 



Type 


Label 


Mass 


Spins 


SO{6) x 50(3) 


Bosonic SO (6) 


X 


6^3 


< j < k - 1 


(6,2j + l) 


Bosonic 50(3) 


a 


3^3 

i 

3 


< j < k - 2 
l<j<k 


(l,2j + l) 
(l,2j + l) 


Fermionic 


X 
V 


4^3 
12 ^ 3 


2 < 3 < k - | 
\ <3<k-\ 


(4,2.7 + 1) 
(4,2j + l) 



Table 1: Oscillator masses for k coincident five-branes (adapted from pf) 



The contribution of the oscillators labeled by u x" is 

k-l 



g 0,2, + D^+s = e^ + f-ta + i^ + t^-i)^) m 

3=0 X ' 



which converges in the large k limit (as x < 1) to 

6(x 2 + x 6 ) _ 6x 2 (l + x 4 ) 



(l-£ 4 ) 2 (1-x 4 ) 2 



(3.4) 



The computations for the other kinds of oscillators are similar, and we will only display the 
large k limit result. The "a" oscillators contribute x 4 (l + x 4 )/(l — x 4 ) 2 , the "/?" oscillators 
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(3x 4 -x 8 )/(l-x 4 ) 2 , the "x" oscillators 8x 5 /{l-x 4 ) 2 and the "77" oscillators 8x 3 /(l-x 4 ) 2 . 
The sum of all those contributions should equal one, so after multiplying by (1 — x 4 ) 2 we 
get 

(x + l) 4 (x 4 -4x 3 + 4x 2 -4x + l) = (3.5) 
having a unique real root between zero and one, 

x = 1 + -J= - ^V2 + ^ ps 0.323556 (3.6) 

so = - log A • ps 13 ^ 406 and T H ps 0.0738519cj. Note how close these values are to 
those of the single five-brane (k = 1) case. 

3.3 Membranes 

The matrix model vacuum in which the SO (3) matrices are in the N dimensional irreducible 
representation corresponds to a spherical membrane. The vacuum with k copies of the n 
dimensional irreducible representation, with n = N/k — > 00, and k fixed, corresponds to k 
coincident spherical membranes || . The spectrum of excitations above those vacua was 
computed in || || and is given in table |2| (again in units of u). 



Type 


Label 


Mass 


Spins 


50(6) x 50(3) 


Bosonic 50(6) 


X 


6^3 


< j < n - 1 


(6,2j + l) 


Bosonic 50(3) 


a 




3^3 
3 
3 


< j < n - 2 
1 < j < n 


(l,2j + l) 
(l,2j + l) 


Fermionic 


X 
V 


4^3 
12 ^ 3 


\<3<n-\ 


(4,2j + l) 
(4,2j + l) 



Table 2: Oscillator masses for k coincident membranes (adapted from B) 



In the single membrane case, the gauge group is completely Higgsed, and the oscillators 
are not matrices anymore. Although there is an infinite number of such oscillators as 
TV — ► 00, their density is just linear (the mass grows as the spin j while the degeneracy is 



proportional to 2j + 1), and as explained in subsection 2.3, the Hagedorn temperature is 
infinite. 

In the case of k coincident membranes, the gauge group is Higgsed down to U(k), and 
the oscillators are k dimensional matrices. Each oscillator, then, is comprised of k 2 entries, 
a fixed finite number, and viewing those as the basic oscillators, even disregarding the 
gauge invariance condition, shows that the density of oscillator masses is still insufficient 
to produce a Hagedorn behaviour. 

4. Thermodynamics of the Matrix Model 

In the previous section, we demonstrated that the spectrum of the matrix model has a 
Hagedorn density of states in the large iV limit and that we should expect a phase transition 
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at some temperature of order the parameter u. In this Section, we will demonstrate 
that this is indeed the case using a direct analysis of the thermodynamics of the matrix 
model. We will be interested in the thermodynamic partition function, where we take the 
Hamiltonian, H, for the model corresponding to Ql.2|) and form the partition function 



Z[N,R,f3] = Tre-P H = e "W^ 

where (3 = 1/T is the inverse temperature and F[N,R,/3] is the free energy. Before we 
discuss this, we shall again illustrate the method by using a simple toy model. 

4.1 The Toy Model 

Consider again the toy model of section 0. The thermodynamic partition function has a 
functional integral representation 

Z = J [dA][dX j ]e-^H^Mi D ^ 2 +^) (41) 

where the fields have periodic boundary conditions in Euclidean time 

A(t + P) = A(t) , Xj(t + ff) = A>(r) (4.2) 

To proceed, we must fix a gauge. The periodic boundary condition prevents us from 
choosing the A = gauge. The best we can do is to fix a gauge where the gauge field is 
static and diagonal, 

A -A(t) = , A ab = 5 ab A a (4.3) 



( DXj ) ab = —Xf + i(A a - A b )Xf (4.4) 



dr 

In that case, 

d_ 

dr 

The Faddeev-Popov determinant for fixing this gauge is computed in the Appendix. It 
turns out to be similar to the Vandermonde determinant which appears in the measure of 
unitary matrix integrals, 

detDdet'4- = T\sm-\A a - A b \ (4.5) 
dr -"-jL 2 

a=f=b 

where det is a determinant restricted to the space of non-constant modes. Then, the 
integral over the oscillator fields in (4.1) is Gaussian and can be performed explicitly 



Z= j l[dA a l[ S m^\A a -A b \l[det-i (- 

J 0=1 a^b ' a,b V 



dr 



A h 



(4.6) 



The determinant can be evaluated explicitly (see the Appendix, also ref. |^]). This model 
with d = 1, and a different critical behavior from the one which we shall be discussing 
in the following, has been studied extensively in the context of lower dimensional string 



theories. For example, see refs. [23|, 
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We obtain 



Z= f J{dA a J[ S iu^\A a - A b \\{ ( 

J a=l a^b " a,b \ 

The effective action for the eigenvalues is then 



sinh ( § (w + i(A a - A b )) 




S e s = Y,\ dln ( sinh f ( u + - A b))) - Insin^ - A b \\ (4.8) 

a^b ^ ^ ' > 

We have dropped an A-independent ground state energy. There are now N remaining 
degrees of freedom, the variables A a , and because there are two iV-fold summations, the 
action is of order iV 2 . For this reason, in the large N limit the integration in ( ]j.6|) can 
be done in the saddle point approximation. In that case there is a well-defined expansion 
of the integral as a power series in 1/N 2 . To begin, we must find the saddle point by 
minimizing the effective action. The equation for the minimum, for each a, is 



V d ; T " ~7 ! 1 x = T «* J (A a - A b ) (4.9) 
{-f cosh /3a; -cos/? (A a - A b ) 2 v b> v ' 

b^a b=f=a 



It is useful to introduce the eigenvalue density, 

N 

A" 



1 N 

P(0) = ^y^0-PA a ) (4.10) 



a=l 



and to assume it is a smooth function in the large N limit. Without loss of generality, 
—tt<8 = (3 A < 7r, and the density is normalized so that 

d0p(6) = 1 (4.11) 

Obviously, the density is constrained to be positive, p{9) > 0. Now, the action can be 
written as 

S cS =N 2 J d9de'p{6) P (9') Idhx (^sinh + 1(0 - 9')j^ - In fsin^|0 - 0'M j (4.12) 

We see there is a sort of competition between the eigenvalue repulsion of the second term, 
coming from the Vandermonde, and the eigenvalue attraction of the first term, coming from 
the action. The attraction is more pronounced in the high temperature regime, /3 — * 0. 
The equation for the eigenvalue density is given by 

dd 'p^ d ,^r} e ~ d 2 m = I <w> cot \ ( e - & ') ( 4 - 13 ) 



cosh /3oj - cos(9 - 9') J ry ' 2 

where it is understood that the principal value around the singularity should be taken in 
the right hand side integral. Note that this equation should be satisfied only in the support 
of p(9), that is, for angles 9 such that p{9) > 0. This equation should be solved together 
with the normalization condition ( 4.11] ) and the positivity constraint. 
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First, we observe that the constant eigenvalue density 



* = h (414) 



is always a solution of ( |4.13| ) where both sides of the equation vanish. However, this solution 
is stable only for a range of values of (3. To see this, we Fourier expand inside ( 4.12; ), 6 



/ °° "I _ r]p-nfiu 
dOdO' p{6)p{0') — cos(n(6 - 6')) (4.15) 
n=l 71 



From this, we can see that if (3 is large enough, then the phase where the density is constant 
(so that the Fourier modes vanish, J d9p(9)e = for n ^ 0) is stable. As the temperature 
is raised, and (3 is decreased, the first instability sets in at [3 = Ph = 1/Th where 

T « = £ <" 6) 

This coincides with the Hagedorn temperature in ( [2.71 ) which we found by estimating the 
high energy density of states in section ^. 

Note that the approximations of section [2] amount to keeping the n = 1 term only 
in ( f4,15| ), and that the arguments there are invalid in the deconfined phase above the 
Hagedorn temperature. However, the result of section [2] correctly gives the free energy 
below the Hagedorn temperature. All that ( 4.15| ) tells us about the low temperature phase 



is that the free energy there is zero in terms of TV 2 ; indeed, it is 0(N°) in the confined 
phase, but more elaborate computations are needed in order to reproduce this from the 
matrix theory perspective. 

The straightforward generalization of this model where we have different frequencies 
for the bosonic oscillators and possibly some fermionic oscillators gives the general equation 
for the critical temperature of the matrix model, 



This also coincides with the equation for the Hagedorn temperature which was found in 
section ||. 

The phase transition that we have identified is the usual Gross- Witten |2f| transition 
of unitary matrix models. In that phase transition the distribution of eigenvalues of the 
unitary matrix rearranges itself from one which has compact support near the identity 
matrix to one which is distributed on the whole unit circle. 

Above the phase transition the eigenvalue distribution is not uniform. Without loss 
of generality, we may assume it is an even function concentrated around 6 = 0. The 
distribution can be approximated by the semi-circle one, 



^)=4 c °4v!- sin2 2 (4i8) 



'Here, we use the identity 



cos n6 



/ . u \ v— v cos nu 
m y sm — j = — h const. 
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This distribution has support in the region | sin || < y | and the parameter £ is restricted 
to the range £ < 2. This approximation is exact when we truncate the low temperature 
expansion of the matrix contribution to the effective action at its first nontrivial, p(9) 
dependent, term 

dln(l - e -P"-W-<>')) „ - de -P" e -K0-(>') (4.19) 
and the saddle point equation is approximately 

J dO'p{e') \ 2de-P" sm{9 - 9') - cot i(0 - 9')^j = (4.20) 
This equation is solved by ( 4.18| ) when 7 



( = 2 1-Jl-l/de-^ (4.21) 



This solution exists only when de~@ w > 1, that is when the temperature is greater than the 
Hagedorn temperature. In this region £ < 2 and the semi-circle distribution is well-defined. 
Note, however, that the expansion parameter in the approximation ( [4.19 ), e - ^, is at least 



d ' 



This approximation is not too bad, however, even in the high temperature limit. For 
[3 = we get a narrow distribution, as £ ~ \- ^ is easy, though, to solve exactly for the 
density in this limit. Equation ( [4.13 ) reduces, for d > 1, to 



/ 



dO'piO') cot \ (9 - 6') = 



Because of the implicit principle value we immediately see that in this limit, p(9) = 5(9) 
is a solution (with support only at 9 = 0); the eigenvalue attraction totally wins over. 

4.2 Confinement-Deconfinement Transition 

The matrix model phase transition that we have found is characteristic of a deconfining 
phase transition in a gauge theory which is heated to a sufficiently high temperature. 



The matrix quantum mechanics has the gauge symmetry (2.2) where, since all variables 
transform under the adjoint representation, the unitary matrices are periodic up to an 
element of the center of the group, 

U(t + (3) = zU{t) (4.22) 



7 To see this, we need the integrals 



dB'p{6') cot -(9 -0') = ^sin( 
2 4 



where again the principal value is implied, and 



J d6p{6) 



cos e = [ i - 4 
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and where z is an element of the center of the group. If the gauge group is U(N), the 
center is U(l). If it is SU(N), the center is Zjv, the multiplicative group of the iV'th roots 
of unity. 

No local operators are sensitive to the center element in ( 4.22j ). However, there is a 



nonlocal operator which transforms under the center, the Polyakov loop operator 

p = Tr V exp (i J dr^(r)^ (4.23) 
where V denotes path ordering. Under a gauge transform with boundary condition ( 4.22; ), 



P -> zP (4.24) 

The center symmetry is often treated as a global symmetry of the finite temperature 
gauge theory. As a global symmetry one can ask the question as to whether or not it is 
spontaneously broken. An order parameter is the expectation value of the Polyakov loop, 
(P). If the center symmetry is unbroken, this expectation value must vanish. If it is broken, 
then (P) can be non-zero. 

The vanishing or non-vanishing of (P) has an interpretation in terms of confinement. 
The expectation value, 

, v f[dA}[dX}e-Io^ 2 Tr((DX)^X^ TVl ^ dT A(r) 

J[dA] [dX]e~ X dT * Tr(( DX ) 2 +" 2x2 ) 

is the ratio of partition functions for the system where one classical quark source is inserted 
to the partition function in the absence of the source. F — Fq is interpreted as the difference 
of the free energies with and without the quark. When the expectation value vanishes, this 
is interpreted as taking an infinite amount of energy to insert the quark. Thus, the phase 
with unbroken center symmetry is confining. When the symmetry is broken, the free 
energy of the extra quark is finite and that phase is deconfined. Some examples of lower 
dimensional gauge theories where the expectation value can be computed appear in ref. 



1 33j , 34, Ba|. In particular, the case of the non-abelian coulomb gas could be written as an 



effective unitary matrix model and it exhibited a phase transition very similar to the one 



that we have studied here [34, p5| . 



4.3 The Matrix Model 

The variables in the action of the matrix model ( |1.2j ) can be rescaled so that it has the 
form 



+ ie abc X a X b X c + ^a a [X a , if>] 
\e a ^ aI gij[X a ' \^ J ] + ±e a/3 r I 9'i'][X a ',^ J ] - \[X\X^A (4.26) 
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Here, a a are usual Pauli matrices and gfj relate the 5t7(4) fundamentals to an S0(6) 
fundamental || . As discussed in the Appendix, rescaling of the coordinates does not affect 
the path integral measure. 

We will consider the semiclassical expansion of this model about the classical vacuum 
corresponding to the single five-brane, 

X cl = 

Expansion to one loop order around this solution gives the semi-classical partition 
function 



z«TT/ [dAg] {WlF{iy 1 \™bw\\™bwu> a (427) 

Jo (det B (-D 2 + 1)) 3/2 (det B (- J D 2 + 1/4)) 3 

Where the differential operators in dets and deti? have periodic and anti-periodic boundary 
conditions, respectively. We discuss how to evaluate the determinants explicitly in the 
Appendix. The result is 

z « TT A"d TT M-ltfr+^-^miiA.-Ai 

Ll J« tfi, |Bin i + - A)) | 3 |sm | (& + i(A a - A t ))\" 

When N is finite, this integral is a smooth function of pi/T. When TV is put to infinity, 
it can have singularities and it typically has a phase transition at some critical value of 
p/T. In that limit, the action, which is the logarithm of the integrand, is of order N 2 , 
whereas there are ./V degrees of freedom, A a . In the infinite N limit, the integral can 
be evaluated by saddle point technique. When /i/T is a large number, the integrand is 
dominated by the Faddeev-Popov term which tends to spread the eigenvalues apart, giving 
a uniform distribution of e lAa around the unit circle in the complex plane. When \ijT is 
small, on the other hand, the terms in the denominator cancel the Faddeev-Popov term 
and the remainder of the integrand emphasizes configurations where the eigenvalues are 
close together. Because of the translation invariance of the potential, which is a result of 
the center symmetry of the action, the accumulation point of the eigenvalues is arbitrary 
and breaks the Zjy symmetry of the effective theory. 

In order to evaluate the integral in the infinite n limit, we must find the minimum of 
the action 



d6p{6)d6'p{6') ~ [l + 8(-l) n e-^ n/4T - 3e"^ /3T - 6 e -^ n/6T 

n=l n 



e in(0-e') 



It is easy to see that, if the temperature is large enough, the potential favors a non-zero 
value of the Fourier transforms 

p n = r d0p{9)e m6 

J —IT 

In fact, the first one to condense is p\ (po = 1 is constrained by normalization). What 
keeps it finite is the constraint that p(9) > for all values of 9. 
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The critical temperature for p n satisfies the equation 

I _|_ g^_l)"g-W 4T crit _ 3 e -M™/3T crit _ Q e -tm/6T clit _ q (4.29) 

At the phase transition the free energy changes from a function which is zero in leading 
order to one which is non-zero. The phase transition is first order. 

The reason why one can have a phase transition in such a low dimensional system (the 
only dimension is the one labeled by 8 which comes from the matrix index) is because the 
effective action is non-local, having infinite range interactions. 



5. Discussion 



5.1 M-Theory Perspective 

The energy scale u of the excitations above the different vacua was recognized in Q to be 
the scale fi appearing in the M-theory plane wave background and in the matrix action 
( |1.2| ). The higher order interactions are governed by the coupling 

where we work in Planck length units, l v = 1. Moreover, in the five-brane case, the energy 
shifts of the excitation masses from the free, quadratic, values are given by the 't Hooft 
coupling 

\ = &N=(j?) N (5.2) 

so our results are valid only when A is small, A < 1. Obviously we need also to take 
N — > oo for the 't Hooft limit in the matrix model. Furthermore, we need to keep /U fixed 
in order to get a nontrivial Hagedorn temperature result. All this means that we take the 
limit in which the DLCQ radius vanishes, R — ► 0. The brane light cone momentum, 

4- N ,~ s 

P R (5 ' 3) 

therefore diverges in this limit, p + —* oo. The radius of the five-sphere on which the five- 
brane wraps can be computed in the small coupling regime [[l(J. Classically the SO (6) 
directions are zero, but quantum mechanically, at one loop order, they give 

-w) w -(^r-(^-T-(^< <m. 

so it diverges if A is kept fixed, but not necessarily if, additionally, A is taken to zero in 
this limit. 

The five-branes we have studied are therefore quite peculiar from the M-theory point 
of view. As solutions with a light-like time direction, they are in the infinite momentum 
frame, moving (almost) at the speed of light. They have large momentum on another 
light-like direction which is compact and small. This makes it difficult to compare them 
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directly with ordinary NS five-branes or M five-branes. The five-branes might be flat (with 
r$ large in Planck length units) or not, according to the exact limit taken. 

We would have liked, of course, to explore the behaviour of the five-branes in the large 
N limit where \i and p + (and therefore r§) are fixed. This however implies R — * oo, so 
(7 e fT — * oo and, a fortiori, A — > oo, where our approximations are invalid. 

Anticipating the results of the next subsection, let us look at the five-sphere radius 
in terms of the length scale ^ (which is the typical (classical) length scale of our matrix 
model, as can be seen from the rescaling done in ( 4.26| )). We obtain 



depends only on the 't Hooft coupling and is directly proportional to it. The five-brane 
radius r 5 c j can also be derived from the the classical five-brane action |1C], 

This calculation is valid for large values of A. measuring again in terms of ^ we get 

This again depends only on A, but has a different power dependence. It is interesting to 
note that for the crossover region A ~ 1, r§ and r 5 c j are of the same order. It is natural 
to assume that there is a smooth transition between the small and large 't Hooft coupling 
results. 

When comparing a single five-brane and a stack of k of those, we would like to deal 
with the same radius, rs, and momentum per brane, p + ; as we have seen, both must be 
large. The equations ( |5.3| , [5~4| ) still apply, provided we replace N by n = N/k in them; n 
is, so to speak, the amount of "iV per brane" . We therefore get that [i is fixed, and hence, 
as we have seen in this section, the Hagedorn temperature is fixed up to a very small 
numerical factor. As we need n fixed as we vary k, the total iV is proportional to k, so we 
are not directly comparing different vacua of the matrix model with the same parameters. 
This result should be contrasted with the Little String Theory case, arising from a stack of 
flat coincident NS five-branes in Type II string theory, where the Hagedorn temperature 
behaves [16, [l?]] as 

In addition, from subsection |2.2.4 we know that there is no polynomial prefactor for 
the exponential Hagedorn density of states. 

For the membrane case, the constraints given by the validity domain of the approx- 
imations are much more lax, and our results are more robust. It is not entirely clear 

3 /2 

H whether the corrections to the free, quadratic, results are governed by (^jjj^ or by 
N 1 / 2 i-^jj) ■ In any case, that magnitude might be small while R being either large or 



."V 

small in the N — > oo limit. The momentum of the membrane, p + , and its (classical) radius 

1 uN , . 

n = ^ (5.8) 
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should be large, although the ambiguity just mentioned does not allow us to decide whether 
they diverge in the large N limit or whether they can be held fixed. 



5.2 String Theory Perspective 

The DLCQ of M-theory, giving matrix models, can be thought of as the limit of quan- 
tizations of M-theory on nearly light-like circles, when an auxiliary parameter R s with a 
dimension of length vanishes, R s — ► 0. By boosting, this parameter can be taken as the 



radius of a spatial compactification of another M-theory 21]. This latter compacti- 
fication gives rise to a type IIA string theory on a certain background. Our aim in this 
subsection is to correlate the small R s limit with the large N limit so as to obtain a valid 
and illuminating string theory description of our system. 

The original M-theory background is a plane wave with parameter //, and it is com- 
pactified on an (almost) light-like circle with radius R. Performing an (almost) infinite 
boost, it can be regarded as an M-theory background of a plane wave with a different 
parameter fj, s , compactified on a spatial circle with radius R s . 

More precisely, we start from a plane wave background 

ds 2 = -2dx' + dx'~ - ((y) 2 (x a f + (^) 2 (x a ') 2 ) {dx' + f + (dx a f + (dx a ' f 
F+123 = lis (5.9) 
with an identification 

' ' J f ' \ 

\ x' 10 + 2ttR s 




(5.10) 

\ x j \ x -\- ah rt s I 
where x /=t = \{t' + x /10 ). 8 This is related to our original background ( |1.1[ ) by the boost 

(t \ | t' cosha — x' 10 sinha \ (5 11) 

x 10 J I — t' sinha + x' cosha I 

or 

x + = e - a x' + ,x~ = e a x'~ (5.12) 
and the plane wave parameters are related as 

ii = e a fi s (5.13) 

The identification in the (t, x 10 ) coordinates is 

t — 2ttR s sinh a \ 
+ 2ttR s cosh a I 

+ 2kR ( ~} ) +0{e~ a ) (5.14) 





8 This convention is diflterent from that of the other sections by a factor of y/2, which is not relevant for 
our rough estimates. 
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where 

R = e a R s (5.15) 

Thus the light-like compactification x~ ~ x~ + 2ttR with a plane wave parameter \x can 
be regarded as a limit of the spatial compactification x' 10 ~ x' 10 + 2-7ri? s in the infinitely 
boosted frame a — * oo. In this infinite boost limit, R s is taken to be zero as R s ~ e~ a so 
that R remains finite. 

From ( 5.13j ) and ( 5.15| ), and reinstituting the Planck energy scale in the original M- 
theory, we obtain 

RM 2 R S M 2 

v - = ^ (5.16) 

/" Ms 



in the a — > oo limit. Instead of changing fi to fi s , we prefer, following [19]-[21], to rescale 
the Planck length. Put otherwise, we map the original M-theory to a new M-theory, 
compactified on a spatial circle with radius R s , having a plane wave parameter /i (not fi s ), 
and a Planck energy scale M p . This new theory is denoted by M, and all of its parameters 
will be written with a tilde. The advantage in doing so is that the light cone energy scale 
of the excitations in the original M-theory (which is [i in our case) is kept fixed under this 
map. Explicitly, we demand 



and therefore obtain that 



R s Ml _ RMp _ R S M% 



M 2 M 2 
p 



(5.17) 



(5.18) 



Ms M 

or, from ( |5.13| ), that M 2 = e a M 2 . To keep M p finite and fixed, we need M p — > oo. We 
also rewrite (|5.17[) as 

R S M 2 = RM 2 (5.19) 

The M theory is described by N DO-branes in weakly coupled type IIA string theory, 
with [21] the string coupling 9 

~g s = Rl'\RM 2 p fl A -» (5.20) 

and string energy scale 

M s = J R; 1 /4( i?M 2 ) 3/4 (521) 
The 't Hooft coupling for the M theory is seen to be the same as the one in the M-theory, 

f3RM 2 \ 3 [3R S M 2 \ 3 f3R s M 2 \ 3 

r=h^ r=h^r (5 - 22) 



and in string units it is expressed as 



X = g sN l^\ (5.23) 



9 These relations for flat space are valid only for a weak background. This will turn out to be the case, 
see below. 
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In the standard prescription of DLCQ in refs. |19|-[ffifl, R s is simply taken to be zero. We 
wish to consider a version of DLCQ which is appropriate for our 't Hooft limit. We scale 
R s in an N dependent way in the large N limit, so that the 't Hooft coupling g s N of 
open string theory on DO-branes is kept finite. From ( 5. 23] ) this means M s j 'fi is kept finite, 
which is different from the case studied in refs. [|l9|-(2l]]. Later we will see, however, that 
although finite the string energy scale should be kept large compared with jj,. Then from 

o 

(M±j „ \V%NR s y 1/4 (5.24) 
so we take the following large N scaling of R s to keep M s //i finite, 

R a „JL(-£-) (5.25) 



This is equivalent to taking the boost parameter as 



R s \^Ml ( a \ 4 



R V 2 \M, 

From ( 5.22]) we also need to take 



iV 2/3 (5.26) 



The radius r§ of the spherical five-brane (|5.4| ) is also rescaled, as a transverse direction, in 
the new theory, 

r b M p = r 5 M p (5.28) 

Then 

fl = -^r 2 »(-£-) 5.29 

Measuring the rescaled radius in the string units, we obtain 

^~A(4-) 2 (5.30) 

Thus the radius of the transverse spherical M five-brane, or equivalently spherical NS 
five-brane in IIA string theory, is kept finite and small in string units. 

We need to check the validity of the flat space formulas ( ^20|) and flUf) we have used. 
The plane wave metric in (t, x 10 ) coordinate is 

ds 2 = ~{l + ^-)dt 2 + ^(1 - ^){dx 10 ) 2 - ^-dtdx 10 + {dx a f + {dx a ') 2 , (5.31) 
where 

2 



F 2 



(I) (xa)2+ © ( * a ' )2 (5 - 32) 
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Thus the flat space relations flpg ), (|5T2l|) are valid when the transverse length scales we 
are considering are much smaller than thus we need 



fJ>r 2 , fir 5 < 1 (5.33) 

where 

r 2 (L- (5.34) 

R S M$ 1 ' 

is the size of the fuzzy membrane solution. This is the case when the 't Hooft coupling is 
finite and 

„ f2 ~(£) 2 «l (5.35) 

i.e. fi/M s is kept finite and small. This condition justifies the use of the DO-brane action 
without stringy correction. 

We conclude that our 't Hooft limit is the limit where the string scale and the radius of 
the spherical NS five-brane are fixed. String theory in an NS five-brane background with 
a fixed string scale gives the little string theory (LST). The validity of our calculations 
is limited to small A where the NS five-brane is highly curved in string units, and this 



prevents us from direct comparison with the results of LST on flat space [18]. However, 
the existence of the phase transition appears robust and is expected to occur over a range 
of A values. 

To recapitulate, the original M-theory compactification has the following parameters: 
the Planck energy scale M p , N, the null circle radius R, and the plane wave parameter \i. 
We look at the compactification as the limit R s — ► of nearly light-like compactifications, 
and map, via a boost, to a string theory with the following parameters: the string energy 
scale M s , N, the string coupling g s and [i. Taking the R s — > limit in a certain way, we are 
able to satisfy the assumptions underlying our approximations: in the M-theory side, those 
are the large N limit, N — > oo, and the weakly coupled matrix model condition A < 1; we 
take the 't Hooft limit of fixed A. In the string theory side, there are the decoupling limit 
condition g s — > and the weak background condition \i -C M s near the origin of the space, 
where the branes reside. We then conclude from ( |5.30| ) that the NS five-brane size is finite 
and small in string units, 

r 5 M s < 1 (5.36) 

5.3 Outlook 

If we do not fix the light-cone momentum, the finite temperature partition function of the 
matrix model should use the rest frame energy 

P° = ^={p + +p-) (5.37) 

in the Boltzmann weight for the partition function, 

Z = Tte ~M p++ P~) (5.38) 
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The trace over p + is a sum over N and the trace over the states of p is the partition 
function of the matrix model, 



Z=J2e-$*e-f ,F V , M (5.39) 

N=0 



Here, F[[3, N] is the free energy of the matrix model. The convergence of this sum for large 
N is determined by the free energy in the large N limit. 

In the membrane states, the free energy is of order one in the large N limit and the 
partition sum converges. In this same limit, the M five-brane states are strongly coupled 
unless we take the 't Hooft limit. In that case, there is a phase transition and the behavior 
of the free energy changes from being of order one to order N 2 . Further, when it is of order 
iV 2 , it is negative and the partition sum diverges. 

This divergence could be interpreted as an instability of the finite temperature system. 
The logical speculation is the collapse of the system to a black hole. This is certainly what 
we would expect in finite temperature field theory on a Minkowski background. There, we 
would expect that the Jeans instability renders any theory of quantum gravity unstable at 
finite temperature. Indeed, in the BFSS matrix model, we do expect the free energy to 
grow faster than linearly in TV. 

A. Fixing the Static Diagonal Gauge 

The thermodynamic partition function is defined by the Euclidean path integral 

Z = J[dA]...exp J^ /T dr^TT{(DX) 2 + ...) ) j (1.1) 

where T is the temperature (we use units where Boltzmann's constant ks = 1) the variables 
have periodic boundary conditions, 

A(t) = A(t + 1/T) , X(t)=X(t + 1/T) , ... 

This model has the gauge invariance, 

X^U(t)X{t)U\t) , A(t)^U(t)(A(t)-i^Jtf(t) 

where U(t) = zII(t + 1/T) is a periodic function of time and z is constant element of the 
center of the group. It is present here since the gauge transform is entirely in the adjoint. 
In a sense, the true gauge group is the factor group U(N)/Z^. There is a quantity, the 
Polyakov loop operator 

P = Tr (Ve 1 ti /T drA ( T )^j 

which transforms non-trivially, 

P -> zP 

This is interpreted as a global symmetry of the theory. Its realization is related to the gauge 
interaction in the theory and can be used to analyze confinement. There are two phases. 
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When the symmetry is good, the Polyakov loop operator averages to zero: (P) = 0. When 
the symmetry is spontaneously broken, the Polyakov loop operator can have non-vanishing 
expectation value, (P) / 0. Of course in this very low dimensional system, a discrete 
symmetry can only be broken in the large N limit. 

The expectation value of P is interpreted as the energy of this system with one addi- 
tional external fundamental representation source. Indeed, the path integral quantization 
of the system with such a source would replace the partition function Ql.ip by (P) and 
with n such sources the partition function would be replaces by (P n )- Thus, we interpret 

F[T] = -T In ((P)) 

as the free energy that it would take to introduce a fundamental representation source 
into the system. When (P) = this free energy is infinite and we say that the system is 
confining, whereas when (P) ^ this symmetry is finite and we say that the system is 
deconfined. 

Let us analyze the partition function (1.1) in more detail. We can fix a static gauge 
by enforcing the gauge condition 



d 
7fr 



A(t) = 



Then, taking into account the Faddeev-Popov ghost determinant, the partition function 
becomes 



Z 



[dA(T)}...8(dA/dT)det 



expU^ 1/T driTr((DX) 2 + ...) 



The measure in the path integral can be defined using the Fourier transform, 

oo 

A(t) = VT A n e 27TinTr 

n=— oo 

The integration measure can be defined in this basis, 

[dA(r)} = HdA n 

n 

The Faddeev-Popov determinant is a determinant over the non-zero modes 
5(dA(r)/dT) det 



IT II 6 ( 27rinTA n) ■ (2mnT)(2mnT + iA^) 

a,b n^O 



= Hl[5(A n ).[(2iTinT + iA^) 

a,b n^O 

where A^ is the adjoint. Integrating the modes of A n with n / and decomposing the 
measure into an integral over the eigenvalues of A and the unitary matrices necessary to 
diagonalize yields the measure 



iV 



[dA] = \\dA a \{{A a - A b ) 2 [dtt] 

b+a 



0=1 
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where is the unitary matrix which diagonalizes A. Also, when A is diagonal A^ 
A a — A\y. Then, the total measure becomes 

N oo / \ N 

n dAa n n ^ inT + ^ - a >^ n ^ inT 

o=l 6^a n=-oo \n^0 

iV 00 

= nylin (27rinr + <(A,-^)) 

a=l by^an=—oo 

where we have used the zeta-function regularization to evaluate the infinite product 

00 

J] 2mnT = W^nT) 2 = {2nT) 2 ^e^'^ = i 

n^O n=l 

Riemann's zeta function is defined by 



00 .. 

«•>-£;? 

n=l 

and has the values 

C(0) = ~2 > C'(0) = --ln(27r) 

To take the remaining product, we use the formula 
00 / 
]7 (2mnT + i<f>) = exp I ^ ln(2vrmT + ■ 

n=— 00 \ n 

= ex P (-|#|^j; (27rmT 1 +#)2 

= exp (- / #/ # J- £ dz± coth(,/2T) I -^ 

where we have taken two derivatives in order to get a convergent series and the contour C 
encircles all of the singularities of the cothz/2T function located at z = 2mnT. Doing the 
contour integral by collapsing the contour onto the singularity at z = — i(f) gives 

f[ (2mnT + i<f>) = exp(-J dcp j d^^—^ coth(-^/2T)^ 

= exp (^J d4>(^cot(<f>/2T)+p)} = exp(lnsin(0/2T) + a + 0<f>) = e a+ ^ sm{4>/2T) 

where a and (3 are arbitrary constants parameterizing the subtractions that are needed to 
renormalize the product. If we arbitrarily set both of them to zero, we get the gauge-fixed 
measure for integrating over gauge fields as 

The absolute value sign is there because the Faddeev-Popov determinants should always 
have absolute value signs. 
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